Driving quantum walk spreading with the coin operator 
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We generalize the discrete quantum walk on the line using a time dependent unitary coin operator. 
We find an analytical relation between the long-time behaviors of the standard deviation and the coin 
operator. Selecting the coin time sequence allows to obtain a variety of predetermined asymptotic 
wave-function spreadings: ballistic, sub-ballistic, diffusive, sub-diffusive and localized. 
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I. INTRODUCTION 

The quantum walk (QW) has been studied as a nat- 
ural generalization of the classical random walk in rela- 
tion with quantum computation and quantum informa- 
tion processing [l], Q ■ The possibility that future quan- 
tum algorithms will be based on the QW has attracted 
the attention of researchers from different fields. This 
topic is a matter of continuous attention, as an exam- 
ple it has been recently suggested by Childs [|[ that QW 
can be regarded as a universal quantum computation. 
Some possible experimental implementations of the QW 
have been proposed by a number of authors 0, H, IS 0] j 
beyond them, the fact remains that the development of 
experimental techniques to trap samples of atoms using 
resonant exchanges of momentum and energy between 
atoms and laser light [8| may provide a realistic frame to 
implement quantum computers @. 

One of the most striking properties of the one- 
dimensional QW is its ability to spread over the line 
linearly in time as characterized by the standard devi- 
ation <j(t) ~ t, while its classical analog spreads out as 
the square root of time a(t) ~ t 1 / 2 . 

Several authors have studied the QW subjected to dif- 
ferent types of coin operators and/or sources of deco- 
herence to analyze and verify the principles of quantum 
theory as well as the passage from the quantum to the 
classical world. In Ref. 10] it was shown that the phases 
in the coin operator and the initial state of the coin can 
be used to control the evolution of the QW. The pres- 
ence of decoherence in the QW has been studied as one 
possible route to classical behavior [III, [I^, HH • On the 
other hand the appearance of a small decoherence can be 
used to enhance some properties of the QW in the devel- 
opment of quantum algorithms [Hj]. The QW subjected 
to multiple independent coins was analyzed in Ref. [l4j, 
showing that this type of excitation can lead to a diffu- 
sive spreading instead of the ballistic one. The QW with 
a position dependent phase is studied in Refs. [H, [l6j] 
where dynamical localization was found. A coin with ex- 
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plicit time-dependence is introduced for the first time in 
Ref. [13], also finding dynamical localization and quasi- 
periodic dynamics. In Ref. [l8| a non-linear dependence 
on the position-chirality probabilities is introduced in the 
evolution, finding a variety of dynamical behaviors, in- 
cluding ballistic motion and dynamical localization. The 
QW was also generalized [I|| introducing two coin op- 
erators arranged in quasi-periodic sequences following a 
Fibonacci prescription, leading to a sub-ballistic wave 
function spreading, as shown by the power-law tail of the 
standard deviation [cr(i) ~ t c with 0.5 < c < 1]. Simulta- 
neously, we studied the QW subjected to measurements 
[20l | and decoherence [2l[ with a Levy waiting-time dis- 
tribution, finding that the system has also a sub-ballistic 
spreading. 



The appearance of this type of spreading shows the 
existence of a wealth of quantum behaviors that are not 
restricted to be diffusive (decoherent) or ballistic (co- 
herent). These are quantum coherent intermediate sit- 
uations that could be useful to control quantum infor- 
mation. Also, to drive the speed of the wave function 
spreading could be a tool for the development of quan- 
tum algorithms. In Refs. [liJUH some analytical results 
were obtained about sub-ballistic behavior in a stochas- 
tic frame, but they cannot be extended to deterministic 
cases such as in Ref. . Then it remains to clarify how 
some deterministic sequences of the coin operator lead to 
a behavior that is neither diffusive nor ballistic. Here we 
develop a simple QW model with a generalized coin, that 
allows an analytical treatment, to clarify the connection 
between the time evolution of the coin and the type of 
spreading of the system. We show how to select the se- 
quences of the coin operator to obtain a predetermined 
power-law distribution. 



The paper is organized as follows. In the next section 
we develop the QW model with a time depended coin, in 
the third section numerical results are presented, in the 
last section we draw the conclusions, and in the appendix 
A we show the analytical deduction for the moments. 
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II. QW WITH TIME DEPENDENT COIN 

A. The standard QW 

The standard QW corresponds to a one-dimensional 
evolution of a quantum system (the walker) in a direction 
which depends on an additional degree of freedom, the 
chirality, with two possible states: "left" \L) or "right" 
\R). The global Hilbert space of the system is the tensor 
product H S ®H C where H s is the Hilbert space associated 
to the motion on the line and H c is the chirality Hilbert 
space. Let us call T_ (T+) the operators in H s that 
move the walker one site to the left (right), and \L)(L\ 
and \R)(R\ the chirality projector operators in H c . We 
consider the unitary transformations 

U(9) = {T-<E>\L)(L\ + T+® \R) (R\}o {I ® K{6)} , (1) 

where K(9) — <j z e~ l0ay , / is the identity operator in H s , 
and <jy and <j z are Pauli matrices acting in H c . The 
unitary operator U(9) evolves the state in one time step 
t as |*(t + t)) = U(9)\^(t)}. The wave vector can be 
expressed as the spinor 



!*(«)> = E 



k= — oo 
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b k (t) 



\k), 



(2) 



the time t in Eq.® and retaining the terms up to the 
first order in r we have 

„dafc(t) . d9 

2 r smfi — aM) — r cost' 

at at 



(ak+i — flfc-i) cos # sin # 

-— r (a fe+ i cos - ansm ( 



(7) 



In order to make analytical calculations we shall restrict 
to a continuous coin function 9(t) that satisfies the con- 
dition dO(t)/dt -> for t -> oo. Then Eq.® reduces 
to 



dak cos 9 

— (Ofc+1 - ak-i) 

at t 



We define now the effective dimensionless time 

t 

t* = \ J cos 9dt, 

to 



(8) 



(9) 



where t > t >> r and to is chosen in such a way that 
starting from it, Eq.© reproduces correctly the asymp- 
totic behavior of the original map Eqs. (|3l4p . Then to is 
taken as the initial time of the process. Changing to the 
variable t* Eq.® becomes 



where the upper (lower) component is associated to the 
left (right) chirality. The unitary evolution implied by 
Eq.® can be written as the map 



ak (t + r) = ak+i(t) cos 9 +bk+i(t) sin( 
bk{t + r) =Ofc_i(t) sin 6* - b k -i(t) cost 



(3) 
(4) 



} dak 
J dt* 



afc+i — ak-i- 



(10) 



Note that this is the recursion relation satisfied by the 
cylindrical Bessel functions, thus the general solution of 
Eq. (|10p , introducing the initial amplitudes a® — ai(ta) = 
ai(t*=0), is 



B. Generalized QW 

Here we are generalizing the QW to the case where 
different quantum coins are applied. In particular we 
consider a deterministic angular time dependence 9 = 
9(t) for the coin operator. 

In order to uncouple the chirality components in 
Eqs. (|3l4p we consider two consecutive time steps and 
then the original map can be rearranged to obtain 

a k (t + t) sin — (t — r) sin 9 
= afc+i(t) cos# sin#_ — ak-i(t) sin#cos#_, (5) 



bk (t + t) sin 9- — bk{t — t) sin 9 
= bk+i(t) cos9 sin^- — bk-i(t) sin0cos6*_, (6) 

where 0_ = &(t — r). Both components of the chirality 
satisfy the same equation, then in the following we shall 
work with cife(t) only. We are interested in a time scale 
much larger than r. Then making a development around 



a k (t* 



(-l) k - l af J k -i(n 
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where J; is the Ith order cylindrical Bessel function. In 
a similar way, the solution for the same equation for b(t) 



b k (t* 



J2 (-i) k ' l bf j k -i(n 



(12) 



l=-c 



with 6° = bi(to). The initial conditions of these equation 
are not necessarily the same as those used in the discrete 
map Eqs.([51 2]), because the approximation of a differ- 
ence by a derivative does not hold for small times. In 
this context a long time implies many applications of the 
discrete map. The results Eqs. pTl IT2")l . with appropriate 
initial conditions, provide a good asymptotic description 
for the dynamics of the discrete Eqs.® |4|), as we shall 
see below. 

Additionally Eqs.fTTJ [T2"|) show that the long-time 
propagation speed of the probability amplitudes is given 
by v(t) = dt*/dt = cos9/t. 
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The probability distribution for the walker's position 
at time t* is given by 



Pk(t*)~ \a k (t*)\ 2 + \b k (t*)\ 2 , 
that can be expressed as 



(13) 



P k (t*) ~ ]T (-l)" W+0 + &°&°*) Jk-l{f)Jk-j{f)- 

(14) 

We calculate the first and second moments of this distri- 
bution using the properties of the Bessel functions (see 
appendix A), obtaining: 

oc 

mi(t*) ~ -t* ® K a °-i + b °j b T-i] + ™i(0), ( 15 ) 



m 2 (t* 



(**) 



* \2 



1+ £ fc^ + ^J 



3=— oo 

+m 2 (0), 



(16) 



where -ft [a;] is the real part of x and mi(0) and 1112(0) 
are the moments at time t* = (t = to). Note that, if 
we take symmetrical initial conditions, the first moment 
vanishes due to the the symmetrical dependence on the 
initial condition (22J. The asymptotic behavior of the 
standard deviation a = yj mi — 



m\ is 



a(t* 



Bt* + C, 



(17) 



where A , B and C depend on the initial conditions. 
From Eq. (fl"7|) we conclude that if the effective time t* 
remains bounded for all t, then the standard deviation 
remains bounded too. In this case the distribution does 
not spread and is therefore localized. On the other hand, 
when t* is unbounded a increases with t* and the sys- 
tem diffuses. Then, in the asymptotic regime, with t* 
unbounded, we substituted the Eq. pT|) by the following 
equation 



a(n) 



At* 



(18) 



C. Choosing the coin operator 



Among the alternatives for the time dependence of 0(t) 
we study the case 



cos0(t) = 



V2 \t 



(19) 



where a > 0. Substituting Eq.(|T9")) into Eq.([H]) and the 
resulting expression into Eq. (|18p for a < 1 and using the 



Eq. p7|) for a > 1, one obtains 



o-(n) 



TT3 3 yv4(« 1 " Q -"o 1 ^),if«<l, 



4ln(-^ , if a= 1, 
VC, if a > 1, 



(20) 



where we have introduced the discrete time dimensionless 
time n through the definition t = (n — l)r and to = 
(no — 1)t. The standard deviation, in the asymptotic 
limit n » no >> 1, is given by 



if a < 1, 



a(n) 



(21) 



Inn, if a = 1, 
1, if a > 1, 

where we have omitted multiplicative constants. This 
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FIG. 1: The standard deviation a(n) as a function of the 
dimensionless time n in log-log scales. The values of a vary 
from (top) to 0.9 (bottom) in steps of 0.1. The curves for 
the ballistic (a = 0) and the diffusive (a — 0.5) cases are 
shown in thick line to help visualization 

result predicts five different types of asymptotic behav- 
iors: 

a) ballistic for a = 0, 

b) sub-ballistic for < a < 0.5, 

c) diffusive a — 0,5 

d) sub-diffusive for 0.5 < a < 1, 

e) localized for a > 1. 

The result for a = is to be expected because then the 
system reduces to the usual QW with a Hadamard coin. 
In the other cases the system shows a variety of behaviors 
where the standard deviation has a slower growth. The 
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FIG. 2: The standard deviation u(n) (in linear scale) as a 
function of the dimensionless time n (in logarithmic scale). 
The upper curve for a = 1 satisfies the law cr(n) ~ Inn for 
large n. The lower curve for a = 2 oscillates around a constant 
value. 

most extreme of them arises for a > 1. In this case the 
system does not spread when t — > oo and it maintains a 
localized position. 

It should be noted that when the derivative of the coin 
function satisfies d9{t)/dt < r for t — > oo, Eq.© still cor- 
rectly reproduces the asymptotic behavior of the original 
map. Indeed, the approximation made to obtain Eq.Q 
remains valid. In particular, if we take the linear time 
dependence — 27t~/t, with 7 irrational and 7 < t/(2w), 
the present model gives similar results as those found in 
Ref.|l7|]. 

III. NUMERICAL CALCULATION 

In order to verify the approximations made in our an- 
alytical treatment we shall compare the result of Eq. I|2ip 
with the numerical evaluation using Eqs.([31 2]). We take 
as initial conditions a walker starting from the central 
position |0) with chirality [1, i] T ; this choice produces 
a symmetrical evolution as in the usual QW. 

We then calculate the standard deviation a (n) using 
the original map. The results for a from to 0.9 are 
shown in Fig.[TJ where for large n (t > 1000 r) the curves 
have no perceptible difference with n 1- ". 

In Fig. [5] the upper curve, that corresponds to the 
standard deviation for a = 1, shows some quick oscil- 
lations governed by an average logarithmic evolution as 
was numerically checked for large times. The lower curve, 
for a = 2, shows that the standard deviation oscillates 
around a constant value. Then in both cases the average 
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FIG. 3: The probability distribution Pk as a function of the 
dimensionless spatial position k at time t — 5000 r. The val- 
ues of the parameter are, from top to bottom, a — 0, 0.3, 1, 
and 2. The scales for k are different in each case 

behavior agrees with the theoretical predictions. 

The profiles of the distribution Pk for several values of 
a are shown in Fig. [3] at the same time t = 5000 r. For 
a = the known QW profile for the Hadamard coin is 
obtained. In this case the function spreads with its great- 
est speed arriving up to the positions ±3500. For larger 
values of a the spreading speed decreases, the distribu- 
tion shrinks and the two extreme peaks of the distribu- 
tion come closer. When a = 0.5 the standard deviation 
spreads out as a(t) ~ i 1 / 2 , but the probability distribu- 
tion Pk is not Gaussian, it has the same characteristics as 
shown in Fig. [3] for a = 0.3. This confirms that the evo- 
lution corresponds to a coherent unitary process and not 
to the classical random walk. In the extreme case when 
a > 1 the distribution is restricted to small values around 
its initial value and the two extreme peaks have melt in 
a very narrow peak. We conclude that in all cases the 
asymptotic behaviors have an excellent agreement with 
the theoretical predictions. 



IV. CONCLUSION 

In the usual QW the standard deviation has a linear 
growth, however using a biased quantum coin arranged 
in aperiodic sequence a sub-ballistic behavior is obtained 



5 



[19}|. Here we generalize the discrete QW using a time 
dependent unitary coin operator. We find an analytical 
expression for the dependence of the standard deviation 
on the coin operator in the asymptotic regime. This ex- 
pression shows explicitly how to select the unitary evolu- 
tion to obtain a predetermined asymptotic behavior for 
the wave function spreading: ballistic, sub-ballistic, dif- 
fusive, sub-diffusive or localized. The parameter a of 
Eq. (fT9)) determines the degree of diffusivity of the sys- 
tem for large times. This feature opens interesting pos- 
sibilities for quantum information processing because it 
could be used for controlling the spreading of the wave 
function. 

Finally another point worth mentioning is that the 
asymptotic differential Eq. (flO|) has the same form that 
the differential equation for the amplitudes of the quan- 
tum kicked rotor (QKR) in resonance (23|. The simili- 
tude between the QW and the QKR in resonance, sug- 
gested in Refs. [2(| HH , has been recently extended to 
the anomalous diffusion behavior presented in this work 
[25l | . Therefore, as the QKR has been realized experi- 
mentally [26| and the quantum resonances have already 
been observed [27], [28| , our results for the QW could be 
experimentally confirmed through an implementation of 
the QKR in resonance. 

We acknowledge the support from PEDECIBA, ANII 
and thank V. Micenmacher for his comments and stimu- 
lating discussions. 



APPENDIX A 

The first and second moments of the distribution Pk 
are 



nil 



n(t) - fcPfe W' ( A1 ) 

k— — 00 
00 

m 2 (t) = J2 k2p k(t)- (A2) 

k— — 00 

Substituting Eq. (fT4|) into these equations we have 

00 

mi(t) = ]T (-l)-M(c!la?+b?b?)F Jl ,(A3) 

j,l=00 
OO 

m 2 (t) = £ (-1)-^ (a?a°* +bfb°/)G jl ,(A4) 

3,1=00 



where 



Fji = J2 k Jk-i(t)J k -j(t), (A5) 

k— — 00 
00 

G jt = e Jk-i(t)J k -j(t). (A6) 



With the change of indexes fx = k—l and v — I — j , these 
equations become 



J2 nMt)j^- v (t) 

f-i— OO 

OO 

+1 M*)Jm-»®, (AT) 

f_L—~ OO 

OC 

fj.— — oc 

00 

+21 i* Mt)Jn-At) 

f-i— OC 

OO 

+l 2 MQJp-uV). (A8) 



Gji 



In the above equations, three different type of sums are 
involved. They can be expressed in terms of the Kro- 
necker delta using the propertied of the Bessel functions 
(Ref/29.], p. 992, Eq. 8.530) 



^ Jn(t)Jii- v (t) = S„o, 
/1— —00 

00 ^ 

——00 
00 

J2 » 2 Mt)J»-„(t) 



(A9) 



- + , (A10) 



(S v -2 + 2S uQ + 5 y2 ) 



+^{5 v -i-5 u i). (All) 



Substituting the above equations into Eqs. (|A7[ IA8[) and 
returning to the original indexes 



Gji 



|) <y«_ 2 + i(2i + i)<y w _ x 



(A12) 



+ / 2 



+i(2l-l)8 lj+ i+(\) 6 h .j. (Al:-!) 



Finally, Eqs.([T5l[TH]) are obtained substituting Eqs. (|A12[ 
IAT3) into Eas. (jA3llA4jl . 



fc — — OO 
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